Zero-bias tunneling anomaly in a clean 2D electron gas caused by smooth density 

variations 
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We show that smooth variations, Sn{r), of the local electron concentration in a clean 2D electron 
gas give rise to a zero-bias anomaly in the tunnel density of states, even in the absence of 

scatterers, and thus, without the Friedel oscillations. The energy width, ujo, of the anomaly scales 
with the magnitude, 5n, and characteristic spatial extent, D, of the fluctuations as {Sn/ D)'^^^ , while 
the relative magnitude Si/fv scales as {5n/ D). With increasing uj, the averaged 5u oscillates with uj. 
We demonstrate that the origin of the anomaly is a weak curving of the classical electron trajectories 
due to the smooth inhomogeneity of the gas. This curving suppresses the corrections to the electron 
self-energy which come from the virtual processes involving two electron-hole pairs. 



PACS numbers: 71.10.Pm, 73.40.Gk, 73.23.Hk 
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Introduction. The origin of a zero-bias anomaly in the 
tunnel density of states of disordered metals had been 
traced [Til to the enhancement of the electron-electron 
interactions, caused by their diffusive motion. In two 
dimensions, the relative correction, 5v{uj)/v, to the tun- 
nel density of states due to this enhancement is equal 
to {l/2iTE^T)\n{E^T'^ /oj)\iy{u}t) [2i]. Here v = m/irh^ is 
the bare density of states, is the Fermi energy, m is 
the electron mass, and r is the scattering time. Diffusive 
description applies in the energy domain u) < 1/t. In 
clean samples with mobility ^10^ cm^/V s this domain 
is very narrow, ^ 10~^meV. In fact, as it was demon- 
strated in Ref. y, the 2D zero-bias anomaly extends into 
the ballistic regime U) ^ 1/r, and essentially retains its 
functional form. Virtual processes, responsible for the 
anomaly in this regime, involve one impurity and one 
electron-electron scattering with either small, q ^ kp, or 
large, q ~ 2kF, momentum transfer. 

The relative magnitude of the interaction correction, 
bvjv, falls off with increasing the electron mobility. As 
experimental samples become progressively cleaner, the 
question arises whether the tunnel density of states in the 
absence of impurities exhibits a zero-bias anomaly. This 
issue was first addressed in Ref. iJ; the calculation in this 
paper predicted the interaction correction of the form 
5v\uj)/v cx UJ. However, later analysis [sj indicated that, 
for a finite interaction range, d, the singular behavior, 
Sv/v = oj/AEy, of the correction saturates at w < Wp/d, 
where is the Fermi velocity. 

In the present paper we identify a new mechanism of 
a zero-bias anomaly, which is at work for finite-range 
interactions and in the absence of impurities. Namely, 
we show that a narrow feature in (5z^(u;) emerges as a 
result of weak, long-scale, variations of the electron den- 
sity, nir), which are generic for high-mobility samples. 
Our main idea is that the high-order electron-electron 
scattering processes in a clean 2D gas, i.e., the processes 
that involve more than one virtual electron-hole pair, are 
anomalously sensitive to the variations of n{v). An ex- 
ample of such process with two virtual pairs is shown in 
Fig. la. The diagram in Fig. la with three interaction 



lines describes creation of an electron-hole pair, which is 
subsequently rescattered into another pair, and, finally, 
annihilated. As was first pointed out in Ref. 0, the mo- 
menta of states, involved in this process, are strongly 
correlated, namely, they are either almost parallel or al- 
most antiparallel to each other. It is this correlation that 
is affected by the spatial inhomogeneity. The resulting 
suppression of the contributions of the higher-order pro- 
cesses, like shown in Fig. la, to the self-energy, gives 
rise to a zero-bias anomaly. Lack of strong correlation 
in the momenta directions in excitation of a single pair 
implies that second-order processes do not contribute to 
the anomaly. 

Qualitative consideration. The degree of alignment of the 
momenta of states in the diagram Fig. la can be estab- 
lished from the following consideration. Denote with Q 
and P the momenta transfer in course of creation and 
subsequent rescattering of the electron-hole pair. Then 
the conditions that the energies of all electrons and holes, 
constituting the pairs, are close to the Fermi surface can 
be presented as |eqj+Q/2 - Ep\ ~ |eqj_Q/2 - E^\ ~ uj, 
and |eqi-p+Q/2 - E^\ ^ |eq^+p_Q/2 - E^\ ^ uj, where 
tq — q^ /2m, and uj is the energy of the pair. The above 
conditions can be met when either |Q| and |P| are both 
small (much smaller than k-p), or when one of them is 
small, while the other is close to 2^^. For concreteness, 
we consider the case |Q| « 2kp, |P| ^ fcp. Then it follows 
from the first condition that |qi | ^ k-p, and that qi • Q ^ 
ujk-p/v-p. Similarly, the second condition requires that 
(P ± qi) • Q ^ ujkp/vp. Combining the two conditions, 



Q/2+q, Q/2+qi-p 




FIG. 1: (a) A diagram describing a virtual process of creation, 
rescattering and annihilation of the electron-hole pair; (b) 
Illustration of lifting the momenta alignment due to curving 
of electron trajectories in external field. 
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we have |eqj_|_p_Q/2 — eqi-Q/2| ^- The latter relation 
can be cast in the form P • (2qi + P — Q) ~ ujkp/v^. 
Since both scalar products P • Q and P • qi are ~ Lukp/v^, 

1 /2 

we arrive to the estimate |P| ~ |qi| ~ k^{uj/E^) ' . 
Therefore, the angle between the momenta within the 

1 /2 

first pair is small as Iqil/fcp ~ (uj/Ep) . Similarly, the 
momenta within the second pair are aligned in the angu- 
lar interval ~ (uj/ E^)^^^ . 

For the purpose of our derivation, we reformulate the 
above restriction in coordinate space, where co defines the 
distance, r, between the subsequent scattering processes 
via the relation uj ~ Vp/r. Correlation between the mo- 
menta implies that the three points, r = ri, r = 0, and 
r = r2, in which creation, rescattering, and annihilation 
take place, are located close to the same straight line, 
see Fig. lb. The "tolerance" in the angle between the 
vectors ri and r2, is the same as the degree of alignment 

1/2 

in the momentum space, 9{r) ^ (l/fcpr) ' . 

In the presence of inhomogeneity, the Fermi momen- 
tum, fcp = (27rn)-'^/^, becomes a function of coordinates. 
It is convenient to characterize the random spatial varia- 
tions of n(r) by a random force field, F(r), related to the 
local density gradient as Vn(r)/(n) = eF{r)/Ep, where 
(n) is the average density. Denote with D ^ k^^ and 
6n <^ (n) the characteristic scale and the magnitude of 
the density fluctuations. Then the typical value of the 
force is F ^ {E^ / eD) (Sn/ (n)) . The force, F(r) curves 
slightly the classical electron trajectories transforming 
them into arcs with curving angle 6^ = eF±r/2Ej^, 
Fig. lb, where F±_ is the component of force perpendic- 
ular to r. Obviously, the process represented by the dia- 
gram in Fig. la, gets suppressed as 6'f exceeds 9{r). The 
condition 6'f = 0{r) defines the characteristic distance 



k-\E,k,/eFf/^ 
and the corresponding energy scale 

Wo = fp/^o ~ Ep{kpD) 5n/{n) 



2/3 



(1) 



(2) 



The latter scale is the energy width of the feature, di'{uj), 
in the tunnel density of states. As seen from Eq. this 
scale is determined by the characteristics of the density 
variations in combination {6n/ D)^^^ . 

The scales tq and uq can be derived qualitatively from 
a different reasoning. The phase acquired by the electron 
upon travelling the distance r, is 0(r) = k^r. Elongation, 
SC, of the trajectory due to curving, results in additional 
phase 6(j){r) = k^SL = kpr{9p — s\n9p)/9,, ~ kpr9r{r)'^, 
where the curving angle, 9jr{r), was determined above. 
Curving becomes important when 5(j){r) ^ tt. This con- 
dition yields the same r = rg as given by Eq. 

In the above consideration we assumed that the force 
does not change within the characteristic distance, rg, 
between the collisions. The corresponding condition. 
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FIG. 2: Dimensionless function, describing the shape 

of the zero-bias anomaly is plotted from Eq. (|19[) versus di- 
mensionless energy z — uj/uq. Inset in the lower-right corner: 
enlarged plot of $(2:) in the domain 3 < z < 5. 



To <C -D, can be cast in the form 

2/3 ,^2^1/2 



D 



Epkp 



kpD \ eF 



{{n)D^ 



,2/3 



< 1. 



(3) 



Eq. ^ requires that the density variations are very 
smooth, D ^ (n)-^/^/ {6n) . The other point to be 
checked upon is whether the language of the smooth 
variations of local density, n(r), that we have used, is 

adequate. Position-dependent n(r) can be introduced if 

1/2 

the statistical fluctua.tion, ((n)Z?^) , is smaller than the 
change, SnD^, of the number of electrons within the cor- 
relation area, D^, due to the smooth fluctuations. It is 
seen from Eq. ([3]) that our main condition rp <C -D is 

1 /2 

stronger than the condition SnD^ ^ , so that 

the reasoning within the language of local density fluctu- 
ations is justified. 

It is also instructive to compare the width, wqi with 
characteristic spatial change of the potential energy of 
the electrons, U — Epdn/{n). As seen from Eq. ([2]) 



ojq/U - {6nD^)-^/^ < 1, 



(4) 



so that the anomaly is much narrower than the variation 
of the chemical potential. Concerning the magnitude, 
Svq ~ Siy{LL!Q), of the anomaly, we will establish that 




SnD^ 



> 



1 



(5) 



((n)L>2)3/2 ((n)i:'2)3/4 

in course of the analytical calculation, to which we turn. 

Green functions. Finding the functional form of Siy(uj) 
amounts, essentially, to evaluation of the diagram Fig. la 
in the coordinate space with account of the random (but 
locally homogeneous) field, F. This field enters into the 
electron Green function 

Gn(0, r) = -j== exp + ik^r + iS^{0, r) \ (6) 



in coordinate-energy space via the additional phase, 

S(t>{0,r) = / k{r)dl - k^r, (7) 



where fc(r) is the wave vector along the classical trajec- 
tory, connecting the points and r. Suppose that r is 
directed along the x-axis. The parabolic trajectory is 




'XT 



y{x) = Fyx{r - x)/AE^, 



(8) 



while dl = dx\Jl + {dy/dxY ~ dxy I + ^{dy /dx) 
This allows to rewrite Eq. ([7]) in the form 

2 



dx { ^ 

dx 



dx 



k{yix)}-k^ . (9) 
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Substituting Eq. jH]) into Eq. jH) and using the relation 

dk/dy = [k^Fyy{x)]/{2E^), we find 

50(0, r) = ~k^Fy/96El = ~k^r0l/24, (10) 

which, within a numerical factor, coincides with the 
above qualitative estimate. Naturally, the a;-component 
of the field also contributes to Scj). However, this contri- 
bution gauges out in the expression for Sv{ll!). 
Density of states. Analytical expression corresponding to 
the diagram Fig. la in the coordinate space, reads 

liV^ P dfl P 
6„{Lu)^lm-^ J — J dvdv^dv2 Go(r,ri) (11) 

X (ri , r2 )nij_j2 (ri , O)!!^^-^ (0, r2 )G' 
where the polarization operator, nj^(r,r'), is defined as 



FIG. 3: Diagrams representing all third order processes with 
aligned momenta of the virtual states. 



where tq = {2^/'^ /k^) {E^k^ / eFyf'^ . It is seen that 
the characteristic scale of distances ri , r2 in the inte- 
gral Eq. (jl4p is indeed equal to ro in accordance with 
qualitative consideration [see Eq. The origin of the 
combinations rir2(r2 ± ri)/rg can be understood from 
Fig. lb. Scattering sequence ri ^ ^ r2 ^ ri 
leads to the accumulation of the field-dependent phase 
2 [(5(/)(ri) -f (5(/)(r2) — (50 (|r2 — ril)]. The above combi- 
nations emerge from this additional phase upon using 
Eq. (fTO|) . Two contributions to the integral Eq. ([M]) cor- 
respond to the locations of the points ri and Y2 on the 
opposite and the same sides from the origin, respectively. 
The shape of the anomaly. The remaining task is to per- 
form the gaussian averaging over the random field, F . 
Since r^^ cx , this averaging can be performed inside 
the integrand of Eq. (jl4p with the help of the identity 



dx e ^ cos{ax^ + P) — Hi{a)cosl3 — H2{a)smP, (15) 
j-j j-j, (r', r), (12) where the functions Hi and H2 are defined as follows 



and V is the dimensionless (multiplied by v) Fourier com- 
ponent of the interaction potential, which we assume to 
be short-range. We are interested in the oscillatory part 
of polarization operator in the presence of the external 
field. Substituting Eq. (O into Eq. p^ . we readily ob- 
tain for this part 



i?i,2(«) = (7r/2)i/y(l + a2)-i/2 ± (1 + a^)-i, (16) 

We present the final result in the form dviui)/]^ = 
A<^{ijjIujq), with 



nn(0,r) = 



27rr2 



2fcpr - 25(l){v) 



exp 



too — E-p 



e(F 



2\l/2 



V2Epkp 



2/3 



E^iiVn) 



2\l/3 



{Anf' {n) 



(17) 



In Eq. (jlip the integration over azimuthal an- 
gles of ri and T2 can be performed analyti- 
cally, using the relation (exp jp (ri + r2));^p,y^j ,,^^2 = 
sin [p (ri ± r2) + 7r/4] Ip(r\r2)^l'^ . Also the integration 
over r can be carried out explicitly with the help of 
the identity / (irGn(ri, r)Gn(r, r2) = 9Gf2(ri, r2)/917. 
Upon performing these integrations, and combining 
rapidly oscillating terms in the integrand of Eq. (fTT|) into 
"slow", oscillating with period fc^T ^, terms, we obtain 

I r2>r 

X sin [■i;~"'^(ci; — 



in agreement with qualitative estimate Eq. ([2]), and with 
constant A defined as 



A 



. 7 1/2 3/2 



8\/27r3/2 



((Vn)^) 



1/2 



; (18) 



the dimensionless function, ^{z), describing the shape of 
the anomaly, is given by 



dridr2 

2EpT:^ ' '2 kl'^ J r2>r, (rir2)3/2 JO 

^ (^2 ± ri 



dn 



$(z) = 



dpidp2 



(14) 



(piP2)3/2 

P2>Pl 



dz' sin 



{z- z'){pi +P2) 



X sm 



rir-2(r2 ± ri)/r^ + 7r/4 T v^^lu + n){r2 ± n) 



|'5+(Pl,P2)+G+(pi,/)2)+S'-(pi,P2)+G_(pi,p2)} 

= $+(z) + $_(z), (19) 
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where the functions 5*+, C+, and C_ are defined as 
S±{pi,P2) = {pi ± P2)^^^ sin j^{z + z'){pi ± P2) 
x|hi(^PiP2(pi±P2)) - V^j, (20) 

C±ipi,P2) = (pi ±P2)^/^C0S y T (2 + z')(Pl ±P2) 



PlP2(Pl T P2) 



(21) 



In definitions of 5'+ and 5_ we had subtracted from the 
function Hi{a) the zero-F value Hi{0) = ^/tt. Integra- 
tion over z' in Eq. (|19p can be carried out analyticaUy. 
The remaining integrals over pi, p2 were evaluated nu- 
merically. The resulting shape of the zero-bias anomaly 
is shown in Fig. 2. The small-z ^ 1 behavior of $(z) is 
81nz, i.e. it diverges logarithmically. The cutoff is cho- 
sen from the condition that $(z) approaches zero at large 
z. Note, that $(2) exhibits a pronounced feature around 
z = 1. The origin of this feature lies in strong oscilla- 
tions of the integrand in Eq. p4|) . The "trace" of these 
oscillations survives after averaging over the magnitude 
of the random field. In fact, the oscillations persist be- 
yond 2; = 3, as seen in the inset in Fig. 2. Also analytical 
inspection of Eq. for 2: 3> 1 yields 



-23/4, 



sm 



(28/3\/32 



,3/4 



■ exp 



{-28/3z| .(22) 



Other third-order processes. Eqs. (|17[) - (|19p were derived 
for a specific process, illustrated in the diagram in Fig. la. 
However, creation, rescattering, and annihilation of a 
pair can follow a different scenario, e.g., rescattering pro- 
cess can involve the initial electron, as illustrated by the 
second and third diagrams in Fig. 3. Important is, that 
the restriction concerning the momenta alignment, lead- 
ing to the zero-bias anomaly, applies to this scenario as 
well. It also applies to aZZ other diagrams in Fig. 3. Note, 
that diagrams in Fig. 3 do not exhaust possible third- 
order processes @. All contributions to of the dia- 
grams in Fig. 3 have the same analytical structure and 
differ only by numerical coefficients, originating from spin 
degeneracy and from closed fermion loops [each bringing 
a factor (—2)]. Collecting these contributions, amounts 
to multiplying the first diagram in Fig. 3 by 1/2. 
Concluding remarks. Higher order processes in a homo- 
geneous electron gas, involving n > 2 electron-hole pairs 



are also subject to the momenta restriction [6|, leading to 
the anomaly in the presence of inhomogeneity. However, 
these processes are suppressed as (wo/£^f)"^^ due to the 
phase-space limitation. 

Note, that in addition to the oscillating term, the po- 
larization operator Eq. (jl3p contains also a slow-varying 
term cx |f^|/r. Evaluating Si'{uj) from Eq. ([TT|) with this 
slow part of H^ yields a non-anomalous correction. The 
same applies to all diagrams in Fig. 3. Moreover, sum- 
mation of RPA subsequence of the diagrams with "slow" 
no(r), a procedure similar to Refs.13,0, does not produce 
any anomaly, if the interactions are short-range [5| . 

In experimental situations, the electrons are supplied 
to the 2D gas by donor impurities, separated from elec- 
trons by a wide spacer. Growth-related technological in- 
homogeneities, like ridges [7|], do not change the aver- 
age electron density, but redistribute electrons over the 
plane and set scales 6n and D. In principle, indi- 
vidual donors themselves create the Friedel oscillations 
of the electron density. However, these oscillations are 
exponentially suppressed due to the large separation of 
donors from the 2D gas. A question might be asked as 
to why in evaluating of diagram in Fig. la we had ne- 
glected violation of the momentum conservation due to 
the impurity scattering. The answer is that condition 
Eq. ^ justifies such a neglecting. This is because at 
distances larger than the Bohr radius, the donor poten- 
tial is screened. Thus, instead of individual donors, the 
electron is scattered by a smooth in-plane potential with 
spatial scale, D. Then the scattering angle does not ex- 
ceed l/{kpD). On the other hand, the relevant curving 
angle, 9^, is (uJo/E^y/^ - [{Sn/n)Y/^{k^D)-^^^ , as fol- 
lows from Eq. ([2]) . Condition Eq. ([3]) guarantees that this 
angle is bigger than l/(kpD). Summarizing, borrowing 
momentum from donors is not efficient, since they are 
distant and screened. 

This conclusion is also supported by comparison 
of the inverse transport scattering time, t^^^ ^ 
Ep{U/Ep)'^{l/kpD), from screened impurities and the 
energy scale luq. As follows from Eq. ^ the ratio 
l/(woTt,.) - {U/Epy/^{l/kpDy/^ is small. The latter 
also implies that the diffusive anomaly of Ref. develops 
at Lo much smaller than loq. 
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